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Abstract 

Higher transcendental function occur frequently in the calculation of Feynman integrals 
in quantum field theory. Their expansion in a small parameter is a non-trivial task. We 
report on a computer program which allows the systematic expansion of certain classes of 
functions. The algorithms are based on the Hopf algebra of nested sums. The program is 
written in C++ and uses the GiNaC library. 



PROGRAM SUMMARY 

Title of program: nestedsums 
Version: 1.0 
Catalogue number: 



Program obtained from: http : / /www . f is .unipr . it/~stefanw/nestedsums 



E-mail: stefanw@fis.unipr.it 
License: GNU Public License 
Computers: all 
Operating system: all 
Program language: C++ 

Memory required to execute: depending on the complexity of the problem, at least 64 MB RAM 
recommended 

Other programs called: GiNaC 0.8.3, a library for symbolic calculations in C++, is required. It 
is available from |http: //www, ginac . de| . 

Optionally, the program can also be used interactively. In this case the program gTybalt, avail- 
able from [http : / /www . f is .unipr . it/~stefanw/gtybalt .html| is needed. 

External files needed: none 

Keywords: Multiple poly logarithms, Feynman integrals 

Nature of the physical problem: Systematic expansion of higher transcendental functions in a 
small parameter. These expansions occur for example in the calculation of loop integrals in 
quantum field theory within dimensional regularization. 

Method of solution: Algebraic manipulations of nested sums. 

Restrictions on complexity of the problem: Usually limited only by the available memory. 
Typical running time: Depending on the complexitiy of the problem, see also sect. || and |[ 
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LONG WRITE-UP 



1 Introduction 

Higher transcendental function, like hypergeometric functions or Appell functions, occur fre- 
quently in the calculation of Feynman integrals in quantum field theory. Usually they appear 
with a small parameter £ in some of their arguments. For example, within dimensional reg- 
ularization [p, the small parameter £ is given by the deviation of Z)-dimensional space-time 
from the four-dimensional value: D = 4 — 2e. However, a solution in the form of transcen- 
dental functions is not yet useful, since possible poles in £ have not been made explicit. What 
is needed is a Laurent expansion in £. In some cases this can be achieved with the help of an 
integral representation [Q]. However, this method is rather tedious and limited to cases where an 
integral representation is known. Recently, algorithms have been developed for the systematic 
expansion of certain classes of transcendental functions [|^]. These algorithms are based on an 
algebra of nested sums. In this paper we report on the implementation of these algorithms into 
a computer library. 

Although these algorithms allow in principle the expansion to any desired order, they grow 
exponentially with the complexity of the problem. Therefore any implementation has to face 
the fact, that at intermediate steps one possibly deals with very large expressions. Commercial 
computer algebra systems are in these circumstances not always as reliable as one would like 
them to be. Computer algebra systems like FORM [Q] or GiNaC [|]], which have been devel- 
oped within the high-energy physics community, seem to be more appropriate. Here we report 
on an implementation based on the GiNaC library. 

GiNaC is a library written in C-i-i-, which allows the symbolic manipulation of expressions 
within the programming language C-i-i-. There are several advantages to this approach: C-i-i- is 
a standardized programming language and supports object-oriented programming. The compi- 
lation of the programs leads to efficiency in the performance. Furthermore, many development 
tools like for the automatic generation of documentation are available. In addition the source 
code is freely available. 

The program library "nestedsums", which we describe here, extends the capabilities of the 
GiNaC library towards the expansion of transcendental functions. This is a domain, which up 
to now is usually not available within commercial computer algebra systems. 

The paper is organized as follows: In the next section the definitions of nested sums and their 
main properties are briefly recalled. Section ^ gives an introduction to the design of the program 
and aims at readers who wish to understand the source code. Section ^ is of a more practical 
nature and describes how to install and use the program. Section ^ deals with checks that we 
have performed and addresses issues like performance. Some simple examples are given in 
section ^ Finally a summary is provided in section 0. 
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2 Theoretical Background 

In this section we shortly review the algorithms for the systematic expansion of transcendental 
functions. They are described in detail in They are based on a particular form of nested 
sums, which we call Z-sums. Z-sums are defined by 



Z{n;mi,...,mk;xi,...,Xk) = V ^•••t^ (1) 



■ ^1 

■mi ■ ■ ■ -mk 

n>il>!2>-->iA>0 1 ^ 



and form a Hopf algebra. If the sums go to Infinity (n = oo) the Z-sums are the multiple poly- 
logarithms of Goncharov [^: 

Z(oo;mi,...,m^;xi,...,Xi) = Li„^,...,„i (xyt, (2) 

Forxi = ... = Xk = I the definition reduces to the Euler-Zagier sums [|^, ||]: 

Z(n;mi,...,m^;l,..., 1) = Z„j,...,„^(n). (3) 

For n = oo and xi = ... = Xk = I the sum is a multiple ^- value [0: 

Z(oo;mi,...,m;t;l,---,l) = C,{mk,...,mi). (4) 

The multiple polylogarithms of Goncharov contain as the notation already suggests as subsets 
the classical polylogarithms Li„(x) [jT^], as well as Nielsen's generalized polylogarithms [jTI]] 

Sn,pix) = Lii,...,i,„+i(l^^^,;c), (5) 

p-i 

the harmonic polylogarithms of Remiddi and Vermaseren [ p^ 

,...,raj- (■^) Li^^^ ,^^( ^1, 1^ , Jc) (6) 

k-l 

and the two-dimensional harmonic polylogarithms introduced recently by Gehrmann and Remiddi 
[|T3p. Euler-Zagier sums occur in the expansion of Gamma functions: 

r(n+£) =r(i+£)r(n) 

X (1 -f£Zi(n - 1) +e^Zn{n - l)+e^Zin{n + ...+e"-^Zn..A{n - 1)) • 

(7) 

The usefulness of the Z-sums lies in the fact, that they interpolate between Goncharov's multi- 
ple polylogarithms and Euler-Zagier sums. In addition, the interpolation is compatible with the 
algebra structure. 

In addition to Z-sums, it is sometimes useful to introduce as well 5-sums. 5-sums are defined 
by 

x'^ X* 

S(n;mi,...,mk;xi,...,Xk) = (8) 
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The 5-sums are closely related to the Z-sums, the difference being the upper summation bound- 
ary for the nested sums: (/ — 1) for Z-sums, / for 5-sums. It is advantageous to introduce both 
Z-sums and 5-sums, since some properties are more naturally expressed in terms of Z-sums 
while others are more naturally expressed in terms of 5-sums. One can easily convert from the 
notation with Z-sums to the notation with 5-sums and vice versa. 



Since the Z-sums form an algebra, any product of two Z-sums with the same upper summa- 
tion bound n can always be reduced to a sum of single Z-sums. For example: 

Zn{n)-Zi{n) = Z2i{n) +Zi2{n) +3 Zni{n). (9) 

The multiplication property is at the core of all our algorithms. In [|3|] we showed that four 
generic types of sums can always be reduced to single Z-sums. These four types are: 
Type A: 

L I L rr^-'rrZ — 7-tZ(? + o- l,mi,...,m/,xi,...,x/), (10) 

Type B: 

x' r{i + ai+bi£) r{i + ak + bke) . , 

L 7^~r^F7^~^ n^-'-Tv^, — — -Z(j + o- l,mi,...,m/,jci,...,x/) 

y"-' r{n - i + a\ + b\z) T{n - i + a'^, + b'^,z) 



X- 



Type C: 



'n-i + c'Y T{n - i + c\ +d[zy" T{n - i + c[, + d'^,z) 

xZ{n — z -fo' — l,mj, ...,mJ/,Xi, ...,xj/), (11) 



i ' (z + c)™ r(/ + ci + Ji £) ■ ■ ■ r(z + ck + dke) 



x5(z + o,mi,...,m/,A:i,...,A:/), (12) 



Type D: 



-E (:)(-»• 



x^ r{i + ai+bie) r{i + ak + bk£) 



{i + c)" r(/ + CI+ diz) ' ' ' r{i + ck + dke) 
x5(« + o,mi, ...,m/,.x:i, ...,x/) 

y"-' r{n-i + a[+b[e) r(n-z + 4, + &^,£) 



X 



n-i + c')""' T{n - i + c\ +d[z)"' T{n - i -f 4, -f J^/E) 
x5(n — i + o\mi,...,mit,x\,...,x'ii). (13) 

Here, all Uj, a'p cj and c'j are integers, c, c', are nonnegative integers and o, o' are integers. 
For sums of type A the upper summation limit n may extend to Infinity. The program library 
implements the algorithms for these four generic types of sums. 
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3 Design of the program 



C++ allows for object-oriented programming. Wherever it was compatible with the required 
efficiency, we tried to write the program according to this paradigm. In C++ objects are rep- 
resented by classes together with methods, which operate on them. All our classes are derived 
from the fundamental class basic, which GiNaC provides. In the following we describe shortly 
the main classes of the program. The complete documentation is inserted as comments in the 
source code and can easily be extracted and converted to latex or html format with tools like 
"doxygen". 



3.1 Letters 

Expressions of the form 

occur quite often in the algorithms. They are represented by the specific class basic_letter 
which has as data members a scale x, a summation index i, a degree m and an offset c. Two 
instances of this class with the same index / can be multiplied, if they have the same offset or if 
both degrees are integers, as for example in 

x' y' {xyY 

(z + c)'"i (? + c)'"2 ~ (z + c)'"i+'"2' 

^' y ^ 1 1 y ^^^^ 

{i + Cl){i + C2) (C2-Cl) (/ + Cl) (C2-Cl) (/ + C2)' 

In the second line the multiplication reduces to partial fractioning. The multiplication operator 
is overloaded and multiplications are performed whenever this is possible. This is illustrated in 
the following code fragment: 

// construct two basic_letters 
ex 11 = basic_letter (x,ml, c, i) ; 
ex 12 = basic_letter (y,m2, c, i) ; 



// multiplication 
ex 1 = 11 * 12 ; 

// 1 is now basic_letter (x*y,ml+m2, c, i) 
If the offset c equals zero, the object 



(16) 



is represented by the derived class letter. We will later see that Z-sums are constructed from 
these objects. Note that letter is derived from basic_letter and not vice versa. This is in 
accordance with the general rule, that a derived class can be substituted everywhere, where the 
base class is required. A further specialization occurs if the scale x is equal to one. The object 

1 

— (17) 

im 

is represented by the class unit_letter, which is derived from letter. Fig. |T] summarizes 
the inheritance relationships for letters. 
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basic 



basic letter 



letter 



unit letter 



Figure 1 : The inheritance diagram for the class basic letter and its derived classes. 



3.2 Nested sums 

Probably the most important class within the program is the class Zsum and its derived classes. 
Fig. ^ shows the inheritance relationships for the derived classes, which correspond to the 
special cases discussed in eq. ^to eq. ^ A Zsum has the data members n and letter_list. 
The former corresponds to the upper summation limit whereas the later is a list of objects of 
the type letter. The actual expressions for the summation indices occuring in the letters 
need not to be known and are therefore replaced by default values. The multiplication operator 
is overloaded and multiplications are immediately performed, whenever this is possible. The 
following code fragment illustrates this: 

ex 11 = unit_letter ( (ex) 1) ; 

ex Zl = Zsum(n, 1st (11) ) ; // Zl is Z_{l}(n) 

ex Zll = Zsum(n, 1st (11, 11) ) ; // Zll is Z_{ll}(n) 

// multiplication 
ex res = Zl * Zll; 

// res is now Z_{21}(n) + Z_{12}(n)+ 3*Z_{lll}(n) 
A similar implementation for 5-sums is provided by the class Ssum. 

3.3 Gamma functions 

In the expansion of transcendental functions. Gamma functions occur always in ratios. The 
class ratio_of_tgamma represents therefore the expression 

r( 1+^728) Tji + ai+bie) 
r{l+bie) {i + a2 + b2e) ' 

The prefactor r(l +Z72e)/r(l -\-b\Z) avoids an unnecessary proliferation of Euler constants 
when this expression is expanded in £. Apart from the obvious data members al, bl, a2. 
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basic 



Z-sums 




multiple poly logs Euler-Zagier sums 




harmonic poly logs multiple zeta values 



Nielsen polylogs 



classical polylogs 



Figure 2: The inheritance diagram for Z-sums and derived classes. 



b2, index and expansion_parameter this class has two integer data members, order and 
f lag_expand_status. The order-variable specifies to which order the object shall eventually 
be expanded. The class is a "smart" class and expands itself into Euler-Zagier sums only if 
f lag_expand_status is set. The expansion is performed according to eq. ^ 
In addition there is a class list_of_tgamma, which is a container for objects of the type 
rat io_of_t gamma and represents a product of terms of the form as in eq. [T^. 



3.4 Algorithms 

The classes transcendental_sum_type_A to transcendental_sum_type_D implement the 
algorithms for the expansion of the sums given in eq. |TD| to eq. [T3|. These classes do the hard 
part of the work and are kept quite general. Since they are kept quite general, they are not 
necessarily particular user-friendly. It is assumed that a user customizes these classes to his 
own needs and writes his own interface. An example is provided by the interface classes in 
the next subsection. We shortly discuss the class transcendental_sum_type_A. The other 
three classes have a similar structure. The class transcendental_sum_type_A represents the 
expression 

r(i+Ji£) r(i+4£) 
r{l+bie)"'r{l+bke) 

^L c^^m-rr^ ^^oV-rf^ Z ; + o - l,mi, ■■■,m/,jci, . (19) 
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Apart from the prefactor r(l + Ji£)/r(l +Z?i£)...r(l +dk£)/r{l +bk£-) this class corresponds 
exactly to eq. [T^. The prefactor avoids that unnecessary terms involving Euler constants appear 
in the expansion. The data members of this class are: n and index, representing the upper 
summation limit and the summation index; letter, lst_of_gammas and subsum, representing 
the expression x'/ (z + c)™, the Gamma functions and the subsum, respectively. These members 
are of type basic_letter, list_of_t gamma and Zsum, respectively. In addition there are the 
data members expansion_parameter, order and f lag_expand_status. 



3.5 Interface classes 

The classes transcendental_f ct_type_A to transcendental_f ct_type_D provide a spe^ 
cific interface to the classes discussed above. They represent: 
Type A: 

r(Ji)...r(^^„) ^ r{i+ai)...r{i+ak) x' 



r(4)...r«,) f^^Y{i + a\)...Y{i + a',_{) i\ 



(20) 



Type B: 



Y{di)...Y{dn) y y Y{i + Ui) . . .Y{i + Uk) Y{j + bi)...Y{j + bi) 

r(4 ).. .r(<o S ;tj, r(/ T0> r( J + ).. .r(j> 

^ Y{i + i^-ci)...Y{i + j + Cm) x\ 4 ^21) 
Y{i + j + c[)...Y{i + j + c'jn jV 



Type C: 



r(Ji)...r(J„) y y Y{i + ai)...Y{i + ak) Y{i + j + ci)...Y{i + j + Cm) x\x{ 
Y{d[)..Tid:,)^^^^Yii + a[)..T{i + a',)Y{i + j + c'^^ il jV 

Type D: 

Y{dl)...Y{dn) y y Y{i + Ul) . ..Y{i + Uk) Y{j + bl)...Y{j + bl) 

Y{d[)..T{d'^,)toPoni + a[)...ni + 4)Y{j + b[)...Y{j + b\) 

Y{i + j + ci)...Y{i + j + Cm) x\x{ 
Y{i + j + c\)...Y{i + j + c'^_,) i\ j\ 

and are modelled on generalizations of hypergeometric functions, the first Appell function, the 
first Kampe de Feriet function and the second Appell function, respectively. Note that the sums 
start at zero and not at one. As an example we discuss the class transcendental_f ct_type_A. 
The constructor for type A is of the form 

transcendental_f ct_type_A (x, a_num, a_denom, d_num, d_denom) ; 

where a_num is a list containing ai to ak- For example for ^ = 4 one has: 

ex a_num = 1st (al, a2, a3, a4) ; 

Each aj has to be of the form aj = k + se, where k is an integer and £ the expansion parameter. 
a_denom, d_num and d_denom are defined similar. This class has a method 

ex set_expansion ( const ex & eps, int order); 

which expands the object to the desired order in £. In addition there is a constructor, which 
expands the object directly. The classes for type B to type D are implemented in complete 
analogy. 
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4 How to Use the Library 



In this section we give indications how to install and use the program library. Compilation of 
the package will build a (shared) library. The user can then write his own programs, using the 
functions provided by the library by linking his executables against the library. 

4.1 Installation 

The program library can be obtained from 

http : / / www . f is . unipr . it / ~ stef anw/ nestedsums 

It requires the GiNaC library version 0.8.3 to be installed. After unpacking, the library for 
nested sums is build by issuing the commands 

. / configure 
make 

make install 

There are various options which can be passed to the configure script, an overview can be ob- 
tained with . /configure — help. 

After installation, the shell script nestedsums-conf ig can be used to determine the compiler 
and linker command line options required to compile and link a program with the nestedsums 
library. For example, nestedsums-conf ig — cppf lags will give the path to the header files 
of the library, whereas nestedsums-conf ig — libs prints out the flags necesarry to link a 
program against the library. 

4.2 Writing programs using the library 

Once the library is build and installed, it is ready to be used. Here is a small example program, 
which defines two Euler-Zagier sums, multiplies them and prints out the result: 

#include <iostream> 

#include "ginac/ginac.h" 

# include "nestedsums /nestedsums . h" 

int main() 
{ 

using namespace GiNaC; 
symbol n ( "n" ) ; 

ex 11 = unit_letter ( (ex) 1) ; 

ex Zl = Zsum(n, 1st (11) ) ; 

ex Zll = Zsum(n, 1st (11, 11) ) ; 

// multiplication 
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ex res = Zl * Zll; 



cout << res << endl; 

} 

After compilation and linking against the GiNaC library and the nestedsums library, one obtains 
an executable, which will print out 

3*EZ (n, 1, 1, 1) +EZ (n, 2, 1) +EZ (n, 1, 2) 

Here, EZ is the output format for Euler-Zagier sums. 

4.3 Documentation 

The complete documentation of the program is inserted as comment lines in the source code. 
The documentation can be extracted from the sources with the help of the documentation system 
"doxygen" [[T^. The program "doxygen" is freely available. Issuing in the top-level build 
directory for the nestedsums library the commands 

doxygen Doxyfile 

will create a directory "reference" with the documentation in html and latex format. 

4.4 Interactive use 

For small problems it is desirable to avoid the editing/compilation cycle and to use the program 
interactively. In addition one would like the output to be in high-quality fonts, like for example 
TeX fonts. gTybalt [|T5Jl is a program which allows one to do just this. It is based on the C-i-i- 
interpreter CINT [ [T6| ] and uses the TeXmacs editor JTtI ] to display high-quality mathematical 
typesetting. It has build-in support for the nestedsums library. To enable this support one has to 
configure gTybalt with the option 

./configure — with-nestedsums 

gTybalt is freely available under the GNU public license from 



http : // www . f is . unipr . it /~ stef anw/ gtybalt . html 



5 Checks and performance 

The library comes with a testsuite, which performs several tests on the implementation of the 
algorithms. Non-trivial tests are checks of known identities, like for example 

1 +eZi{n)+e^Zn{n) + ... +e^Zii...i(n)] \l -eSi{n) +e^Sn{n) - ... + (-l)V5ii...i(n) 

= l + 0{e^+^). (24) 

This relation provides a check on the implementation of the conversion between Z- and S- 
sums and on the multiplication of nested sums. A second class of checks is obtained for finite 
sums by writing out the sums explicitly and to compare the result of a symbolic manipulation 
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with the result obtained from the explicit expression. Furthermore, identities obtained from 
partial integration also yield useful relations which can be used for checks. In the class of 
hypergeo metric functions, Appell function and Kampe de Feriet functions we have for example 
the relations: 

(c — 1)2^^1(^2,^— l,c — l,x) — (c— 1)2^1 (a, Z?, c — l,x) +ax2Fi{a + l,b,c,x) = 0, 

(Ci - l)3F2{a,bi - l,Z72,Ci - 1,C2,X) - (Ci - l)3F2(«,^l,^2,Cl - l,C2,x) 

+xa — 3F2{a + l,bi,b2 + l,ci,C2 + l,Jc) = 0, 

C2 

(c- l)Fi(fl- l,Z?i,Z?2,c- l,Jci,JC2) - {c - \)Fi{a,bub2,c - \ ,xuX2) 

+xibiFi{a,bi + l,b2,c,xi,X2) +X2b2Fi{a,bi,b2 + l,c,xi,X2) = 0, 
(c- l)5i(ai - \,a2,buc-l,cuxux2) - (c - l)5i(ai,a2,^i,c - l,ci,;ci,;c2) 

+xia2—Si{ai,a2 + l,bi + l,c,ci + l,xuX2) +X2«25'i(ai,a2 + l,&i,c,ci,.x:i,X2) = 0, 

(ci - l)5i(ai,a2,^i - l,c,ci - l,xuX2) - (ci - l)5i(ai,a2,^i,c,ci - l,xi,X2) 
CI2 

+x\ai — Si{a\ + l,a2 + l,&i,c+ l,ci,jci,.x:2) =0, 
c 

(ci - l)F2{a,b\ - 1,^72, Ci - l,C2,Xi,X2) - (Ci - l)F2(a,Z7i,Z72,Ci - l,C2,Xi,X2) 

-\-x\aF2{a+\,bi,b2-,ci,C2,xi,X2) = 0. (25) 

For a given choice of the parameters, each term is expanded in £ independently and one verifies 
that each relation yields zero up to the calculated order. 

In addition we considered the expression 

^ _ r(2m - 2e - vi235)r(l + V1235 - 2m + 2e)r(2m - 2e - V2345)r(l + V2345 - 2m + 2e) 

r(vi)r(v2)r(v3)r(v4)r(v5)r(3m - 3e - V12345) 

^^ r(m-e-V5)r(m-e-V23) 2«-2e-v,2345 f f i4 

r(2m-2£-V235) ^ '''^ hh'^-'^- 
T{ii + V3)r(/2 + V2)r(/i + i2 -2m + 2£ + Vi2345)r(/i +12-171 + £ + V235) 



-^1 



X 



r(/i + 1 - 2m + 2£+Vi235)r(/2 + 1 - 2m + 2£+V2345)r(/i +/2+V23) 

2ra-2e-Vi235 

r(;i + 2m - 2£ - Vi25)r(?2 + V2)r(/i + 12 + V4)r(?i + ?2 + m - e - Vl) 

r(zi + 1 +2m-2£-Vi235)r(z2 + l -2m + 2£+V2345)r(zi +Z2 + 2m-2£-Vi5) 

2m-2e-V2345 



X2 



r(zi + V3)r(z2 + 2zn - 2£ - V345)r(zi + Z2 + Vl)r(zi + Z2 + W - £ - V4) 
r(zi + 1 -2m + 2£+Vi235)r(z2 + l +2m-2£-V2345)r(zi +Z2 + 2z«-2£-V45) 



2m-2£-Vi235_^2m-2£-V2345 r(zi + 2m - 2£ - Vi25)r(z2 + 2m - 2£ - V345) 

"^2 r(zi + l+2m-2£-Vi235)r(z2 + l+2m-2£-V2345) 

r(zi +Z2 + 2m-2£-V235)r(zi +Z2 + 3m-3£-Vi2345 
r(zi + Z2 + 4zn - 4£ - V 12345 - V5) 



(26) 



This expression corresponds to a two-loop integral occuring in the calculation of the two-loop 
amplitude for e+e^ 3 jets. For the specific set of parameters (vi, V2, ...,V5) = (1,1,1,1,1) 
and ( 1 , 1 , 1 , 1 , 2) it has been calculated by Gehrmann and Remiddi [|13]] to order £'^ and order £^ , 
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k 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Form 


< 1 


< 1 


< 1 


< 1 


< 1 


2 


11 


58 


323 


1816 


C++ 


< 1 


< 1 


< 1 


< 1 


< 1 


1 


5 


30 


180 


1066 



Table 1 : CPU time in seconds for the multiplication of two harmonic sums on a 300 MHz 
Pentium II with 64 MB RAM. 



respectively. Our program allows us to obtain the result for arbitrary values for the Vj. For the 
two specific cases we confirm the results of Gehrmann and Remiddi. 

Although the program has been written carefully and many bugs have been eliminated dur- 
ing the debugging phase, there is no guarantee that the program is bug-free. It should be clear, 
that if the program is used to obtain results which are published in scientific journals, it is still 
the responsibility of the user (and not of the author of the program) to make sure that these 
results are correct. 



We would like to give some indications on the efficiency of our algorithms. The inner parts 
of the algorithms are based on the multiplication of nested sums and in this part efficiency is 
mandatory. For the specific case of multiplication of harmonic sums there exists already a pro- 
gram called "summer" which is written in FORM [JTS]]. The comparison between the "summer" 



program and our program for the multiplication of two harmonic sums 

'^ll...lW-'^ll...lW' (27) 

k k 

where each factor has k "one"'s in the index field is listed in table ^. As can be seen from 
table H, our program is almost a factor two faster. However, benchmarks tests should always be 
taken with a certain care. Our program is fast, as long as all intermediate expressions fit into 
the available RAM. No particular memory management facilities have been implemented. If 
intermediate expressions exceed the available memory, the operating system resorts to swapping 
pages, which slows down the execution time significantly. 



6 Examples 

In this section we give two example programs for the expansion of a hypergeometric function 
and a Appell function, together with some indications on the CPU time required for the expan- 
sion up to a specific order. 

The following small program 

#include <iostream> 

#include "ginac/ginac.h" 

# include "nestedsums/nestedsums . h" 

int mainO 
{ 

using namespace GiNaC; 
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order 


1 23456 7 8 9 10 11 


2Fi{a£,b£; 1 — c£;x) 


<1 < 1 1 2 4 8 16 31 59 121 285 



Table 2: CPU time in seconds for the expansion of a hypergeometric function on a 300 MHz 
Pentium II with 64 MB RAM. 



symbol a ("a"), b("b"), c("c"), x("x"), eps("eps"); 
int order = 4; 

// expands 2F1 (a eps, b eps, 1 - c eps, x) in eps 

ex F21 = transcendental_fct_type_A (x, 1st (a*eps, b*eps) , 1st (l-c*eps) , 
1st (l-c*eps) ,lst(a*eps,b*eps) , 
eps, order, expand_status : : expansion_required) ; 

// some polishing 

F21 = convert_Zsums_to_standard_f orm (F2 1 ) ; 
cout << F21 << endl; 

} 

expands the hypergeometric function 2F1 (a£, be, 1 — c£,x) in £ up to a given order. Running this 
program will print out 

1+ (eps^3*b*a*c+eps^3*b^2*a+eps^3*b*a^2) *S (1, 2,x) +eps^3*b*a*Li (3,x) *c 
+eps^2*b*a*Li (2,x) 

which agrees with the known expansion 

2Fi(a£,Z7£;l -c£;jc) = l+ab'Li2{x)e^ + ab{c'Li2{x) + {a + b + c) Si,2 W) £^ + 0(£'^). 

(28) 

The function convert_Zsums_to_standard_f orm brings an expression involving Z-sums into 
a standard form. It first removes all Z-sums with non-positive indices. An example would be 

00 

LioW = £^' = T^. (29) 
/=i ^ ^ 

These sums are rather trivial and can be eliminated. It further collects the coefficients of each 
Z-sum and brings the coefficients into a normal form. Table |^ shows the dependency of the 
CPU time for the expansion of the hypergeometric function on the chosen order. 

As a second example we give a short program, which expands the second Appell function 
^2(1, 1,£;1 +£, 1 -e;x,y) in£: 

#include <iostream> 

#include "ginac/ginac.h" 

# include "nestedsums/nestedsums . h" 
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order 


1 2 3 4 5 6 


F2(l, 1, £;!+£, 1 -£;x,y) 


< 1 < 1 2 9 49 279 



Table 3: CPU time in seconds for the expansion of the second Appell function on a 300 MHz 
Pentium 11 with 64 MB RAM. 



int mainO 
{ 

using namespace GiNaC; 

symbol x("x"), yC'y"), eps("eps"); 

int order = 2; 

// expands F2 ( 1, 1, eps; 1+eps, 1-eps; x, y) in eps 

ex F2 = transcendental_fct_type_D (x,y, 1st (1) , 1st (1+eps) , 1st (eps) , 1st (1-eps) , 
1st (1 ), 1st 0, 1st (1+eps, 1-eps) ,lst(l,l,eps) , 
eps, order, expand_status : : expansion_required) ; 

// some polishing 

F2 = convert_Zsums_to_standard_form (F2) ; 
cout << F2 << endl; 

} 

Running this program will print out 

-(-1+x) " (-l)+2*Li (l,x) * (-1+X) " (-1) *eps-Li (l,y+x) * (-1+x) " (-1) *eps 
which corresponds to the expansion 

1 1 1 

F2(l,l,£;l+e,l-£;x,3;) = + - {2\n{l-x) -x-y))z + 0{z^). 

1 —X 1 —X 

(30) 

Table |^ shows the dependency of the CPU time for the expansion of the second Appell function 
on the chosen order. The algorithm for the expansion of the Appell function is more complex 
compared to the one for the hypergeometric function, a fact which is also reflected in the nec- 
essary CPU time. 

7 Summary 

In this paper we have described the program library "nestedsums". This library can be used 
for the symbolic expansion of a certain class of transcendental functions and can be useful to 
scientists in all fields. In particular these sort of expansions are required in the calculation of 
higher order corrections to scattering processes in high energy physics. The library is written in 
C++ and uses the GiNaC library. 
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